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Abstract 

The B —> decay rate is known at the next-to-next-to-leading order in QCD. It is pro¬ 
portional to aerniuY S'licl has a ±4% scale uncertainty before including the O /ml)) 

electromagnetic corrections. We evaluate these corrections and conhrm the earlier Endings of 
Bobeth et ai. Furthermore, we complete the calculation of logarithmically enhanced elec¬ 
tromagnetic efiects by including also QED corrections to the matrix elements of four-fermion 
operators. Such corrections contain a collinear logarithm ln(m^/m|) that survives integration 
over the low dilepton invariant mass region 1 GeV^ < < 6 GeV^ and enhances the inte¬ 

grated decay rate in this domain. For the low-m^^ integrated branching ratio in the muonic case, 
we find B{B —> XsXt~) = (1-59 ± 0.11) x 10“®, where the error includes the parametric and 
perturbative uncertainties only. For B{B Xse^e~), in the current BaBar and Belle setups, 
the logarithm of the lepton mass gets replaced by angular cut parameters and the integrated 
branching ratio for the electrons is expected to be close to that for the muons. 



1 Introduction 


The inclusive decay B —> with / = e or /x is known to be a sensitive probe of new physics 

at the electroweak scale. Its branching ratio has been recently measured by both Belle [1] and 
BaBar [2], In the low dilepton invariant mass region, 1 GeV^ < < 6 GeV^, the experimental 

results read 

B{B XJ+r) = (1.493 ± 0.504+°;^^J) x 10“® (Belle) , 

BiB^Xsi^r) = (1.8 ± 0.7 ± 0.5) X 10“® (BaBar). 

This leads to a world average 

B{B^XJ+r) = (1.60 ± 0.51) X 10“® . (3) 

Measurements for lower and higher values of mjf are available, too. However, for higher mj^, 
non-perturbative effects of the J/T, T' and higher resonances are sizeable, and the theoretical 
predictions have larger uncertainties. On the other hand, for mj^ < 1 GeV^, the branching 
ratio is determined largely by the contribution from almost real intermediate photons, and 
it contains essentially the same information on new physics as is already known from the 
B —>■ Xs'j measurements. Throughout this paper, we restrict ourselves to the dilepton mass 
region G [1,6] GeV^. 

The experimental errors in the branching ratio are expected to be substantially reduced in 
the near future. On the theoretical side, the predictions are quite well under control because 
the inclusive hadronic B —> Xsi^i~ decay rate for low dilepton mass is well approximated 
by the perturbatively calculable partonic b —> decay rate. Thanks to the recent 

(practically) complete calculation [3-9] of the Next-to-Next-to-Leading Order (NNLO) QGD 
corrections, the perturbative uncertainties are now below 10%. 

The branching ratio is proportional to aem(h) whose scale dependence cannot be neglected. 
Indeed, at the leading order in QED, B{B —>■ Xst^i.~) changes from 1.54 ■ 10“® to 1.65 ■ 
10“® when the renormalization scale of CTem is changed from /i = 0{mb) to /i = 0{Mw)- 
This uncertainty is removed by calculating those QED corrections that are enhanced by large 
logarithms ln(M^/M|), where Mh ~ Mw, rrit and ~ rrib, mu. 

In Ref. [9], the QED corrections to the Wilson coefficients were calculated, thereby giving 
most of the electromagnetic corrections that are enhanced by ln(M^/M|). As a result, the 
authors hnd a branching ratio of 1.56 ■ 10“®, B which incidentally corresponds to setting = 
«L(/^ mb) at the leading order in QED. We have calculated and conhrm the results of Ref. [9] 
for all the two-loop anomalous dimension matrices that determine the size of the ln(M^/M|)- 
enhanced electromagnetic corrections. 

However, there are additional QED corrections that get enhanced by large logarithms, 
namely ln(m^/m^). These corrections are the new result of the present paper. They origi¬ 
nate from these parts of the QED bremsstrahlung corrections where the photon is collinear 

* The number quoted by the authors of Ref. [9] and on which we agree is 1.57 • 10“®. In the text we give the 
result obtained using the updated experimental inputs summarized in Table [T] 
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with one of the outgoing leptons. They disappear after integration over the whole available 
phase space but survive and remain numerically important when is restricted to the region 
that we consider. 

Such logarithmic corrections are found under the assumption that no collinear photons are 
included in the dehnition of the dilepton invariant mass. This turns out to be a very good 
approximation for the muons in the current BaBar and Belle setups [10]. In this case, the 
enhancement of the \ow-m% integrated branching ratio by the collinear logarithms amounts to 
around 2%. The corresponding effect for the electrons would reach around 5%. However, in 
that case, the logarithm of the electron mass gets replaced by the BaBar and Belle angular 
cut parameters and the integrated branching ratio for the electrons is expected to be close to 
that for the muons. We shall describe this issue in more detail in Section El In the preceeding 
sections, our analytical and numerical results will correspond to the case of perfect separation 
of electrons and energetic collinear photons. 

Before we come to the results and details of the calculation, some comments on its systematics 
are in order. Due to the different scales involved, the perturbative corrections come not only 
with increasing powers of some coupling constant, but also with increasing powers of the large 
logarithm L = ln(M^/M|). The perturbative calculation results in an expansion in the product 
of the coupling with L rather than in the coupling alone. 

Because is relatively large, all powers of must be resummed at a given order 

of the perturbative expansion, which is achieved using the renormalization group technology. 
Within this framework, all the logarithms L are absorbed into c* = 0{1). Consequently, each 
electromagnetic logarithm aemL = Csaemlo.s of the conventional perturbative expansion gets 
replaced by f{cs)aern/c(s, where the function /(c^) is found by solving the renormalization group 
equations. Such a replacement of the electromagnetic logarithm is not a matter of convenience 
but an unavoidable consequence of resumming the QCD logarithms and not resumming the 
QED ones. Resummation of the QED logarithms would be technically more difficult and also 
unnecessary, because aemL -C 1. Thus, the conventional expansion in as and Oem is replaced 
by an expansion in ag and m. k = aem/ots- Each order of this expansion is calculated exactly 
in Cg. 

The amplitude oi B ^ Xst^t~ is proportional to ctem- The Leading Order (LO) con¬ 
tributions come from loops and are of order k (the electromagnetic logarithm comes from a 
loop). Higher order terms that are proportional to nag and na"^ are conventionally called the 
NLO and NNLO QCD contributions, respectively. However, since nag = aem, the NLO con¬ 
tributions contain purely electroweak terms, too. Since these NLO terms are enhanced by 
sin^ 6w) while the LO terms are accidentally suppressed, the two contributions turn 
out to be very similar in size. An analogous effect occurs at order the terms of order n'^a] 
are larger than the ones. For this reason, also high terms in the fi:"’a™-expansion remain 
numerically important. 

The corrections to be considered here (and also in Ref [9]) are of order k‘^ and K^ag in the 
decay amplitude. Contributions corresponding to K^a^ ~ in the amplitude will be included 
only if they are enhanced by \n{ml/mj) or by an additional factor of m^/(M^ sin^ 6*iy)- 

Our article is organized as follows. In Section [2] we summarize the results for the branching 
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ratio and explain details of the ^"'a^-expansion. The effective theory used for resummation of 
large QCD logarithms is introduced in Section [H] which is quite technical. It includes the list of 
the relevant operators, the matching conditions for the Wilson coefficients, the renormalization 
group equations and the Wilson coefficients at the low scale. Sections 0] and [5] contain a detailed 
description of the four-fermion operator matrix element calculation. In Section [6] we discuss the 
role of the angular cuts. Master formulae for the branching ratio are summarized in Section [71 
Appendix [^ contains the loop functions that appear in the text. Some intermediate-step 
quantities for the evolution of Wilson coefficients are collected in Appendix [Bl Appendix O is 
devoted to describing techniques that we have used to calculate the QED matrix elements of 
quark-lepton operators. 


2 Branching ratio and numerical results 


In order to facilitate the reading of this rather technical paper, we give the hnal results Erst. 
The differential (with respect to s = ■m.^/'m-b^poie) decay width of i? —>• can be expressed 

as follows: 


dT{B Xsi+j-) 
ds 







(4) 


where the dimensionless function is assumed to include both the perturbative and non- 

perturbative contributions. 

In order to minimize the uncertainty stemming from ml and the CKM angles, we normal¬ 
ize the rare decay rate to the measured semileptonic one. Furthermore, to avoid introduction 
of spurious uncertainties due to the perturbative b XcCu phase-space factor, we follow the 
analyses of Refs. [9,11] where 


C 


Vub ' r(R ^ X,ev) 

Wi, r(R ^ X^eV) ’ 


(5) 


was used instead. Consequently, our expression for the B Xgi'^i branching ratio reads 

2 4 4>u{s) 

C ’ 

where = 1 -|- 0{as, aem, A^/m^) is dehned by 


db[B —>• Xst t ) 
ds 


= B{B^X^e9) 


exp 




Vrb 


( 6 ) 


T{B ^ X^ev) = 


G'im, 


5 

6,pole 


1927r3 


IK/ 4>, 


(7) 


Our expressions for the ratio ^u{s)/^u are summarized in Section [71 Both the perturbative 
and non-perturbative corrections to this ratio are much better behaved than for and 

separately. The factor C = 0.58 ± 0.01 has been recently determined from a global analysis of 
the semileptonic data [12]. All the input parameters that we use in the numerical calculation 
are summarized in Table [H 
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as{M^) = 0.1182 ± 0.0027 [13] 

me = 0.51099892 MeV 

ae(M^) = 1/127.918 

= 105.658369 MeV 

= sin^ 6w = 0.2312 

mr = 1.77699 GeV 

= 0.967 ±0.009 [14] 

me{me) = (1.224 ± 0.017 ± 0.054) GeV [15] 

BR{B Xeez/)exp = 0.1061 ± 0.0017 [16] 

ml^ = (4.68 ± 0.03) GeV [12] 

Mz = 91.1876 GeV 

mt,pole = (172.7 ± 2.9) GeV [17] 

Mw = 80.426 GeV 

mB = 5.2794 GeV 

X 2 — ^ — m%) ~ 0.12 GeV^ 

G = 0.58 ±0.01 [12] 


Table 1: Numerical inputs that we use in the phenomenological analysis. Unless explicitly 
specihed, they are taken from PDG 2004 [18]. 


It should be stressed that the pole mass of the b quark that is present in the dehnition of 
s and in several loop functions gets analytically converted to the so-called IS'-mass before any 
numerical evaluation of the branching ratio is performed. This way one avoids dealing with the 
renormalon ambiguities in the dehnition of the pole mass [19]. The formula that relates the 
pole mass to the IS'-mass can be found e.g. in section 4 of Ref. [20]. 

Let us explain the details of the as and k expansion that we adopt for calcnlating onr hnal 
nnmerical resnlts. The b st^l~ decay amplitnde has the following structure (up to an overall 
factor of Gp)\ 


A = 




Alo + Oig Ajvlo + CKg ^NNLo + G(a^) 

2 yiem ^ 0(a^J 


Al^ + a, 


NLO 


+ ai A 


em 
NNLO 


+ 0{k^ 


( 8 ) 


As mentioned in the introdnction, Alo ~ Anlo and AJ2) ~ C(s Af/^Q. All these terms are 
inclnded in onr calculation in a complete manner, together with the appropriate bremsstrahlnng 
corrections. As far as Annlo is concerned, we use the practically complete results of Refs. [3-9]; 
the only missing parts originate from the nnknown two-loop matrix elements of the QCD- 
pengnin operators whose Wilson coefficients are very small. 

Among the contributions to AJ^j^^q, we inclnde only the terms which are either enhanced 
by an additional factor of mf /sin^ 9w) (with respect to A^^^q) [9] contribnte to the 
ln(mf/m|)-enhanced terms at the decay width level. The latter terms are calculated for the 
hrst time here. They are taken into account in a practically complete manner; the only missing 
part is proportional to the same tiny Wilson coefficient that is responsible for the smallness of 
Alo- 

The pertnrbative expansion of the ratio ^u{s)/^u has a similar strnctnre to that of the 
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squared amplitude: 


= 


K 


+ 


'^lo + “^AloAnLO + i-^^LO + ‘^'^LoAnNLo) 
+a^ 2 {AnloAnnlo + •••) + 

2AlO'^T6 + 2(-^ArLO“^LO + “^LO“^WXo) 

+al 2{Ajyf];^o-^‘]^Lo + -^nnlo-^To + •^lo'^I^nlo) 


+af 2{A^io^7fNLo + ‘^nnlo^^iTlo + 

0{k^) . 


+ 0{at) 


( 9 ) 


In our numerical calculation of $££(s)/$„, we include all the terms that are written explicitly 
in the above equations. The dots at orders and stand for terms that are proportional 
to Alo and A^o and, consequently, can safely be neglected. In the numerical analysis we 
also include subleading l/rric and l/rrib corrections [21-23] as well as hnite bremsstrahlung 
effects [5]. 

Our results for the branching ratios integrated in the range 1 GeV^ < rrij^ < 6 GeV^ read 


Bee 


1.59 ± O.OSscaie ± 0.06,,, ± 0.024c,± 0.015, 


mb 


±0.02c^^(Mz) ± 0.015ckm ± 0.026br,„ 


X 10 


-6 


1.64 ± O.OSscale ± 0.06^, ± 0.025c,rn, ± 0.015, 


mb 


±0.02q,^(Mz) ± 0.015ckm ± 0.026br^, 


X 10 


-6 


(1.59 ±0.11) X 10“® , 
(1.64 ±0.11) X 10“® . 


( 10 ) 

( 11 ) 


The central values are obtained for the matching scale /io = 120 GeV and the low-energy scale 
fib = 5 GeV. The uncertainty from missing higher order perturbative corrections have been 
estimated by increasing and decreasing the scales p.o,fe by factors of 2. Uncertainties induced 
by rrit, rUb, rric, C, as{Mz), the GKM angles and the semileptonic rate are obtained by varying 
the various inputs within the errors given in Table [H we assume the errors on G and me 
to be fully correlated. The total error is obtained by adding the individual uncertainties in 
quadrature. The electron and muon channels receive different contributions because of the 
ln(m^/m^) present in the bremsstrahlung corrections. The difference gets reduced when the 
BaBar and Belle angular cuts are included (see Sec. E])- 

We stress that the indicated uncertainties are only the parametric and perturbative ones. 
No additional uncertainty for the unknown subleading non-perturbative corrections has been 
included. In particular, we believe that the uncalculated order as{fib) ^^^° non-perturbative 
corrections imply an additional uncertainty of around ~ 5% in the above formula. This issue 
deserves an independent study. 

One should also keep in mind that all the effects of the intermediate iji and 'ip' contributions 
are assumed to be subtracted on the experimental side. This refers, in particular, to the decays 
ifi —^ X£~^£~ where low-mass dilepton pairs can be produced. All such decays of the ifi with 
branching ratios down to 10“^ may be relevant. To our knowledge, only X = 7 has been 
considered so far in the experimental analyses. 
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NLO (Q^e?n(/^o)) 

1.81 X 10-6 

NLO (Q^em(/^ 6 )) 

1.68 X 10-6 

NNLO («e„^(/xo)) 

1.65 X 10-6 

NNLO {a,M) 

1.54 X 10-6 

QED (only WC’s) 

1.56 X 10-6 



QED (muons) 

1.59 X 10-6 

QED (electrons) 

1.64 X 10-6 


Table 2: Anatomy of QCD and QED corrections. 


The overall uncertainties in Eqs. CB and fllll) are somewhat smaller than in Eq. (27) of 
Ref. [9]. This is mainly due to the improved experimental value of rrit as well as to our use of 
ml^ rather than mb,po\e- The latter opportunity was already suggested in Ref. [9]. 

In Table [2], we show the partial results that we obtain by adding sequentially all the known 
QCD and QED corrections. The rows denoted by “NLO” and “NNLO” refer to the leading 
order in QED. The row “QED (only WC’s)” contains only those QED corrections that stem 
from the Wilson coefficients. The row “QED” includes all the electromagnetic corrections (that 
are different for electrons and muons, as in Eqs. (ITOil and (ITTll ). 

A numerical formula that gives the branching ratio for non-SM values of the high-scale 
Wilson coefficients of the operators P 7 , Pg, Pg and Pig (see Section [3]) reads 




2.1913 - 0.001655 J(Pio) + 0.0005 I{RwRl) + 0.0535 liRj) + 0.02266 liRrRl) 
+0.00496 XiRjRl) + 0.00513 X{Rg) + 0.0261 X{RgRl) - 0.0118 X{Rg) 

-0.5426 n{Rio) + 0.0281 n{Rj) + 0.0153 ^(PtP^o) + 0.06859 n{RjRl) 

-0.8554 n{RjRl) - 0.00866 n{Rg) + 0.00185 n{RgRlg) - 0.0981 n{RgRl) 
+2.7008 n{Rg) - 0.10705 n{RgRlg) + 10.7687 |Pio|^ + 0.2889 \R^\^ 


+0.00381 iPspT 1-4892 jPgl 


X 10 


-7 


( 12 ) 


a. = 


2.3278 - 0.001655 X(Pio) + 0.0005 J(Pio Rl) + 0.0524 X{Rj) + 0.02266 X{Rj Rl) 
+0.00496 X{R-jRl) + 0.00504 X [Rg] + 0.0261 X{RgR;) - 0.00651 X{Rg) 

-0.5426 P(Pio) - 0.02578 n{Rr) + 0.0153 ^(^7 P^o) + 0-0674 n{RjRl) 

-0.86996 P(P 7 Pg) - 0.0128 n{Rg) + 0.00185 n{Rg P^g) - 0.09926 n{RgRl) 
+2.841 n (Pg) - 0.10705 P(PgPio) + 11.0367 |Pio|^ + 0.2813 IP 7 P 


+0.003765 iPgp + 1.528 |Pg| 


X 10 


-7 


(13) 


where (see Eqs. (fT^ . flUUl) . (170]) for the Wilson coefficient definitions) 


-^ 7,8 — 


^(OO)eff 

^ 7,8 


(+0 


^(00)eff,SM 

'-'7,8 


(+0 


and 


-^9,10 — 


(+0) 

'-'9,10 


(+0 


(14) 
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3 The effective theory 

3.1 Operator basis 

Resummation of large QCD logarithms is most conveniently performed in the framework of 
a low-energy effective theory [24], There are ten operators that need to be considered at the 
leading order in the electroweak interactions. They can be chosen as follows: 

Pi = (sL7/.l'“cL)(cL7^T“6i), 

P 2 = (sl7;,cl)(cl7 ^6 l), 

^3= {sL7^bL)Eq{qi^q), 

Pi= {sa^T%L)Eg{q^'"T'^q), 

(15) 

Ps = J§;^mb{sL(T^^’'T%R)Gl^, 

P9= {sL^f,bL)El{h^l), 

PlO= {sL'^^jbL)Y.l{h^'lbl)- 

In P3, ..., Pe, the quark flavors are q = u,d,s,c,b. In Pg and Pio, all the three lepton flavors 
are present. Contrary to other analyses [25,26], we have not included any gauge couplings in 
the normalization of Pg and Pig. Including them would give only a minor simplihcation in the 
present investigation. 

Once QED corrections are considered, hve more operators need to be taken into account. 
They can be chosen as 

P3Q = {sLif^bL) Eg Qgiqi^q), 

PiQ = {sLl^J^bL) Eg Qg{qYT'"q), 

P5Q= (sL7w7/.2 7/.3^i)Eg<3g(g7^'7'''7^"g), (16) 

PqQ = (SL7/.i7/.2 7M3^“^i) Eg Qgiqi^^Y^Y^T^q): 

Pb = J2 [(sL7Mi7/.2 7M3E)(fe7'''7''"7^'&) - 4 (sl7^6l)(67^6) . 

where Qg are the electric charges of the corresponding quarks (| or — |). 

The Lagrangian of the effective theory reads 


4G r ^ 

Peff — PqcDxQED (u, d, a, c, b, e,M, r) + + Y. C,q{p)P,q + Ci(p)PA . (17) 

V 2 [j=l j=3 



3.2 Matching conditions 

The Wilson coefficients at the matching scale /xq ~ are expanded as follows 

CkifJ'o) = + as(/^o) °^(/xo) + 5s(/xo)^ Cj. °^(/xo) 

+5s(/io)K(/xo) + 5s(/io)^K(/io) + 0{a^, (18) 

where 5^ = a<j/47r. Note, that at the low scale fib ~ ^b, mu, also terms of order k, and H?as 
arise and are included wherever necessary. 

The values of the Wilson coefficients are found from the requirement that all the effective 
theory Green function^ match to the full SM ones at the leading order in At the 

order we consider, the following non-vanishing contributions to Eq. flTSl) must be taken into 
account for the four-fermion operators (s^ = sin^^v^^): 


crGo) 

= 1, 

(19) 



(20) 

o 

o 


(21) 


= I 

3 3 M| ’ 

(22) 

CfVo) 

= jA 1V(I.) - 2Y(x,)] , 

(23) 

Cf>(w) 

= - 2r(x,)], 

(24) 


= -^Y{xt)+ W{xt)+ , 

Sw 9 9 mf 

(25) 

C*io \f^o) 

= — \-y{xt), 

Sw 

(26) 

^3Q^(h'0) 

- 3sf,-9 + 9'”mr 

(27) 

dgAo) 

2 [X{xt) + Y{xt)] , 

(28) 


- 2s2 ’ 

2iSw 

(29) 

Cf°Go) 

= for i = 1, ...,6 , 

(30) 

C'“Vo) 

= Cf‘^\iJ,o) - Cf‘^\iJ,o) for t = 9, 10 , 

(31) 


t For the on-shell IPR functions, the operators from Section ITT] are sufficient. However, it is often more 
convenient find the Wilson coefficients by matching the off-shell IPI functions. Then, additional operators are 
necessary — see Eq. (73) of Ref. [3]. 


9 











cfVo) 


(32) 


x: 


324- 


cifVo) = 


x: 




( 4 s^ - 1) [3 + Tl^\xht) - ^t{Xht) 


3 + Tl,^\xht) - ^t{Xht) 


(33) 


All the one-loop coefficients above have been evaluated in the MS schemeji] The 

functions E{x), X{x), Y{.x)^ W{x), S{x) are collected in Appendix lAl The one-loop coefficient 
is from Ref. [27]. The other one-loop ones have been known since many years (see, e.g., 
Ref. [28]). For and Cf^\ the relevant top and charm (C)^*'^^) contributions to the 

two-loop matching conditions can be found in Section 2 of Ref. [3]. The functions and A*, 
where Xt = /M-^Y and Xht = {Mh/mY^Y, can be found in Ref. [29]. We include also 

the contributions to that were calculated in Refs. [3,30]. Transforming the results of 

Ref. [30] to our operator basis is non-trivialj§ 


3.3 Renormalization Group Equations 

In the effective theory, the RGE for the gauge couplings read 


= -2S‘ Y: Kn,S:ST, 


dog 

dfj, 

dctp 


n,m=0 


“ +2tte X! ^ 


dfi 


n,,m=0 


(34) 


where CKe = The solution for as{fi) with the initial condition at /i = /xo is found 

perturbatively in as(ho) and 5e(/io) but exactly in Vg = 1-|-2/3QQ5s(/io) hi ^ and = 
1 —2/dQQ5e(/io) Including all the 3-loop contributions, and, in addition, the 4-loop pure- 

QCD term, one obtains 


(h) 


«s(ho) _ «s(ho)^ /Y 
Vg vl VA 


s 

10 

•s 

00 


In Vg 


'k 

’00 


hlUe + 


(ho) 


'00 


(1 - Vg) 


+ ( ^ j (in^ Vg -\nVg + Vg- 1) 


A 


+ 


21 

00 , 


In^ Ve + (—2 In Us InUg + pVe InUe) 


YooYi. 


'TO 


+ 


as(ho)^ 


Y30 1 Vs Y20Y10 


A 


00 


(Aoo)^ 


((2n^ - 3) Inn^ + 


+ 



/ 5 

f - In^ Vg + - \iY Vg + 2(1 - Vg) In Us 




^ Beyond tree-level, the Wilson coefficients usually depend on the choice of evanescent operators. Our choice 
is the same as in Refs. [7-9]. 

^ We thank Ulrich Haisch for providing us with the relevant transformation matrices. 
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+ 


®s(/^o) Cle(/^o) 
V^Ve 


Po2 I 1 ^ I 1 '^e Poif^lQ ci I T 'j 

-^[Ve — 1 + ^^PVe In-h , ^ ,n InVe — Wg + 1 


(/3So)‘^ /JSoft'o 


png In 


Inn. 


+ Og X (!?(ttg , Og, OsOe) ) 


( 35 ) 


where p = /3Qoas(po)/(/5oo'^s(po) + /5oo'^e(po))- The corresponding solution for ae(p) can be 
found by obvious replacements: Vg Ve, 5s <->■ 5e and jS-j <->■ —jSfj (also inside the ratio p). 

The MS values of the pure-QCD coefficients /3fg can be found in Refs. [31,32], After sub¬ 
stituting Ca = N = 3, Cp = |) = I cind rif = 5, one hnds Pqq = y, /3fo = P20 — 

and P^Q = ^^|f^C(3) ~ The remaining beta-function coefficients that enter into Eq. 0351) 

reacO 

0§t = = -f. ft'o = I + 3Q?) = f. 

^(;i = (4Cf- 8C^)W = -^ = 4 (O^iV + 3<3?) = f!, (36) 

;3& = f + 2«f^ = If ft'i = 4CFiVg^ = if, 

where Qi = -1, Qu = \, Qd = and = 2Q^ + 3Q^. 


The RGE for the Wilson coefficients reads 

h^<5(p) = 7 ^(p)(7(p), (37) 

where the Anomalous Dimension Matrix (ADM) has the following expansion: 

l{P= E (38) 

n,m=0 

n-\-m>l 

In Eq. fl35|) . we have made no use of the fact that Sg <C Sg. Now we shall take this relation 
into account, and solve the RGE flHTl) perturbatively in 


^ ^ PooMpo) 
~ Poo 5 s (po) 


and a; = 2p^gas(po), 


(39) 


neglecting terms of order Cl(a;^, A^,a;^A^). Let us introduce the following short-hand notation: 


bi = 
64 = 


/^lO 

2(P^o)^’ 

/^ii 

4(/5oo)^/^oo 


bo = 




20 


4(/5oo)'' 


-bl 


bo. = 


A 


01 


2PooPi 


00 


25i53, 


bn — 


Pi 


01 


2 P 00 PI 


-bi, 





{2p^on2PSor 


(40) 


^ All of them except /3fi can be found by modifying the color and charge factors in the pure QCD results. 
As far as f3fi is concerned, we have found it by performing an explicit three-loop calculation. To our knowledge, 
no result for this coefficient has been published before. 
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The known evolution of the gauge couplings 0351) allows us to rewrite the RGE fl37p in terms 
of the variable rj = ag^fio)/ 


^C = - 
drj T] 


ly(io) + ^ ^ik)^k ^ X^u;hW^°^^7]lnr] + C>(o;^ A^ 


k=-2 


c. 


where the matrices are r^-independent 


(41) 


(42) 


B(-^) = CJ + a;2A ,(43) 

(44) 


= ujX{1-, 

= X(i-+ ujX^ + , 


B^'^) = 

The solution to Eq. m reads 


(45) 

(46) 


Cip) = V 


DM+f: ^“’(irt + i: G'">(>)) 

k=-2 k,l=-2 

2 

+ ^ +R{ri)+0{uj\x\uj^x^) 

k,l,m=—2 


where V is the matrix that diagonalizes ITdo) 




V 




E-^G(/io), (47) 


(48) 


The eigenvalues a* and entries of the matrix V are given numerically in appendix [Bl The 
matrices D{t]), and R{r]) depend on the ai and on products E^^'> = 

V-^B^FV . They read 

Ai(h) = Hj\n) = G\f\r]) = Y, E^p^E^pjO^iiiv), 

p 

k7’”’(>))=i:4‘’Gi5d"‘’'"Sf(’)). R.i(r,) = \^ujh{v-'W>°''’v) (49) 

P,1 

The functions fiXiv), dpi'Mt “““d rjf (tj) are given by 

/«Vf) = 


iij o/yi nipj yijy ••'ipqj 

pjai 

1 

aj-\-k— 


- f 

hk-ai yl 


~\~k _ 


when ttj + k — tti = 0, 
otherwise, 


rij\v) = 


^rj°-i In^ 7], 


when ttj + k — tti = 0, 


1 


aA-\-k— 


— _ fiP{r])] , otherwise, 


(50) 

(51) 
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aiplhv) = 


rip\v) 


CL A 


when ttj + I — ttp = 0, 
otherwise, 


In^ T], when ap + k — ai = Ug + I — Up = aj + m — Ug = 0, 


(52) 


h. 


(klm) 

m3 


iv) = 


^ In^r^-r 

l-ap V *<? ^ 

( {k,l+m). . _ ( 

aj+m-ag \^WJ ^ H 


ap+k—o 

_]_ 

aq-\-l—u,p 

1 


(fc) 

ip 


ipq \'U j 5 

(ki) 
ipq 


when Op + fc — Oj 7 ^ 0 and 

Qjg -\- I - Clip = Clj -|- Tfl - CIg = 0, 

when ttq + I — ttp ^ 0 and 

ttj + m — ttg = 0 , 

when ttj + m — ttg ^ 0. 


(53) 


3.4 Anomalous dimension matrices 

In the present Section, we give the ADM’s for the four-fermion operators. When the 
tors are ordered as in the list {Pi,..., Pg, Pg, Pio, -PsQ, ••• , Peg, Pb}, then the matrices that enter 
Eq. (l38ll have the following generic structure: 


opera- 



■ / (nm) \ 

(7CC )2x2 

/ ^ (nm) \ 
[iCP )2x4 

(ffT> 

)2x2 

..(nm) 

17cq 

)2x4 


O 2 XI 



04x2 

(7pp )4x4 

(iff” 

)4x2 

(iff’ 

)4x4 


O 4 XI 


^(nm) _ 

02x2 

[IlP )2x4 

(iff” 

)2x2 

(iff’ 

)2x4 


O 2 XI 

(54) 


04x2 

\1qP )4x4 

(iff” 

)4x2 

..(nm) 

17qq 

)4x4 


O 4 XI 



0ix2 

\1bP )lx4 

(iff” 

)lx2 

..(nm) 

\1bq 

)lx4 

. . (nm) \ 

kIbb )lxl J 


However, the pure-QCD ADM’s have 

a much simpler structure 




( n (nO) \ 

[ICC )2x2 

( 11 (nO) \ 

KlCP )2x4 

02x2 

02x4 


O 2 XI 




04x2 

\1PP )4x4 

04x2 

04x4 


O 4 XI 



^(^0) ^ 

02x2 

02x4 

02x2 

02x4 


O 2 XI 



(55) 


04x2 

\1qP )4x4 

0lx2 (7 

(n0)N 

QQ 14x4 

O 4 XI 




0ix2 

\1bP )ix4 

0lx2 

0lx4 

. . (nO) N 

KIbb )ixi . 



Moreover, four additional blocks vanish in 

7“’ = 0 . iff = 0 . ffp' = 0 . 

We need to know all the non-vanishing blocks of and 


ff”? = 0 . 


( 56 ) 


fflO) 

Tcc = 


■ -4 

8 ■ 
3 

-(10) 

■ 0 

0 

9 

0 ■ 

_ 12 

0 

5 Icp — 

0 

- 0 

3 ^ 

0 _ 


fflO) 

Tpp = 


40 

9 


0 256 


256 

9 


52 

3 

100 

9 

256 

3 

9 


¥ 0 
4 
9 

0 


40 

9 


5 

6 

20 

_2 

3 
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(57) 



-1 

0 

1 

^ico 

0 

0 

_1 


■ 0 -20 0 2 ■ 

0 f 0 0 

-(10) 

40 52 4 5 

9 3 9 6 

0 -if 0 0 

’ "(QQ 

0 -128 0 20 

_ 0 ^ 0 0 _ 


256 160 40 2 

L 9 3 9 3 J 


il,'? = [4], 


-( 10 ) 
Tbp = 


I 0 0 


^( 20 ) 

Tcc 


355 

502 - 


r 1412 

1369 

134 

35 - 

9 

27 

7(20) _ 

243 

243 

243 

162 

35 

28 

) Icp — 

416 

1280 

56 

35 

3 

3 - 


81 

81 

81 

27 - 


^( 20 ) 

IPP 


4468 

31469 

400 

3373 ■ 


■ 832 

4000 

112 

70 

81 

81 

81 

108 


243 

243 

243 

81 

8158 

59399 

269 

12899 

^(20) 

) 7qp — 

3376 

6344 

280 

55 

243 

243 

486 

648 

729 

729 

729 

486 

251680 

128648 

23836 

6106 

2272 

72088 

688 

1240 

81 

81 

81 

27 


243 

243 

243 

81 

58640 

26348 

14324 

2551 


45424 

84236 

3880 

1220 

243 

243 

243 

162 J 


L 729 

729 

729 

243 



404 

9 

2698 

81 

19072 

9 

32288 

81 


3077 

9 

8035 

27 

14096 

9 

15976 

27 


9 

49 

162 

1708 

9 

6692 

81 


1031 

36 

4493 

216 

1622 

9 

2437 

54 


^( 20 ) 

Tbp 


1576 446 172 40 
81 27 81 27 


^,( 20 ) _ r 325 

iBB - 9 


(58) 


Almost all the blocks of are necessary: 


-( 01 ) 

Tbl 


16 

9 


0 


1 


-( 01 ) 

Tcc 



0 

<, 

1- 

r -ii 01 

27 ^ 

.(01) 

tolco 

<i|to 

0 

0 

0 

0 

<1 

'8 -4 ■ 

■0 

-1 

1 

CDlOO ! 

0 

) 7cq “ 

I 0 0 0 

) 111 — 

0 

1 

_ 1 


76 

9 

0 

2 

3 

0 ■ 


16 

9 

0 ■ 


272 

27 

0 ■ 

32 

20 

n 

2 


32 

n 


32 

n 

27 

3 

W 

3 

7(01) 

27 

u 

7(01) 

81 

u 

496 

9 

0 

20 

3 

0 

> 7pl = 

112 

9 

0 

) 7ql — 

2768 

27 

0 

512 

128 

n 

20 


512 

n 


512 

n 

27 

3 

W 

3 


. 27 

u 


81 

u 


(59) 


From 7 ^°^) we need only the mixing of Pi, 


Pe into P9 and Piq: 


■ 11680 416 

-(02) ^ 2187 81 

'CL 2920 104 

729 27 


^( 02 ) 

Tpl 


39752 

136 

729 

27 

1024 

448 

2187 

81 

381344 

15616 

729 

27 

24832 

7936 

2187 

81 


( 60 ) 


The necessary entries of 7 *^^^^ read: 


-( 11 ) 

Tcc 


■ 169 

100 ■ 


' 0 

254 

0 

0 ‘ 


■ 2272 

122 

0 

49 ■ 

9 

27 

-(11) 

729 

-(11) 

729 

81 

81 

50 

8 

) Icp — 

0 

1076 

0 

0 

) ICQ — 

1952 

748 

0 

82 

. 3 

3 . 


243 


L 243 

27 

27 . 


14 





■ 0 

11116 

0 

14 n 


r 23488 

6280 

112 

538 -| 


243 

3 


243 

27 

9 

27 


280 

18763 

28 

35 


31568 

9481 

92 

1012 


27 

729 

27 

18 


729 

81 

27 

81 

7pp “ 

0 

111136 

0 

140 

) 7pq — 

233920 

68848 

1120 

5056 


243 

3 


243 

27 

9 

27 


2944 

193312 

280 

175 


352352 

116680 

752 

10147 


L 27 

729 

27 

9 J 


729 

81 

27 

81 


.( 11 ) 

Tpl 


r 6752 

243 

O' 


r 24352 

729 

O' 

2192 

0 


54608 

0 

729 

^(ib 

2187 

84032 

243 

0 

) 7ql ~ 

227008 

729 

0 

37856 

729 

0 


551648 

L 2187 

0_ 


.( 11 ) 

Tcl 


2272 

729 

1952 

243 


0 

0 


.( 11 ) 

Tbl 


.( 11 ) 

111 



0 16 
16 0 


Finally, the relevant 3-loop anomalous dimensions yield [8] 


.(30) 

Tcc 


12773 

18 

1177 _ 
2 


+ 


1472C(3) 


2144C(3) 


745 _ 4288C(3) 
9 9 


306 - 224C(3) 


.( 21 ) 

Tcl 


' 1359190 I 6976 C(3) q 

19683 243 ^ 

229696 _ 3584 C(3) ^ 

L 6561 81 ^ 


(61) 


(62) 


r 63187 

981796 

202663 

24973 1 


r 1360 

776 

124 

100 1 

13122 

6561 

52488 

69984 

+ C(3) 

81 

81 

81 

27 

110477 

133529 

42929 

354319 

2720 

2768 

248 

110 

L 2187 

8748 

8748 

11664 J 


L 27 

27 

27 

9 J 


(63) 


- (30} _ 

7fp — 


+C(3) 


3572528 

2187 

1651004 

6561 

147978032 

2187 

136797922 

6561 


608 

27 

88720 

81 

87040 

27 

721408 

81 


58158773 

8748 

155405353 

52488 

168491372 

2187 

72614473 


552601 

4374 

1174159 

52488 

11213042 

2187 

9288181 



13122 

6561 

61424 

496 

2821 

27 

27 

9 

54272 

9274 

3100 

81 

81 

27 

324416 

13984 

31420 

27 

27 

9 

166432 

95032 

7552 

81 

81 

27 


6989171 ■ 

11664 
10278809 

34992 

17850329 

2916 

16664027 
17496 . 


.( 21 ) 

Tpl 


(64) 


1290092 I 3200 C(3) 
6561 81 

819971 _ 19936 C(3) 
19683 243 

16821944 I 30464 C(3) 
6561 81 

17787368 _ 286720 C(3) 
19683 243 


0 

0 

0 

0 


(65) 


The three-loop ADM’s have no influence on the logarithmically-enhanced QED corrections at 
the considered order but are necessary for the NNLO QCD corrections. As far as the one- 
and two-loop ADM’s are concerned, we have calculated all of them, and our results agree with 
Ref. [9]. 


3.5 Wilson coefficients at the low scale 

From the solution to the RGE in Section 13.31 we obtain the Wilson coefficients at the scale 
/ift ~ rrih as truncated series in 5s(ho) and We then use Eq. (IHHj) to express the couplings 

at the high scale in terms of as(h6) and For 5^, the simple relation 

5s(/io) = rj 5s(/ife) (66) 
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holds to all orders. In order to obtain the rnnning of k, we invert Eq. 0351) . treating Vg and rj 
as quantities of order 0 ( 1 ), which gives 


^(ho) —-1- 


00 X , Inh 

K {^b) + 


Poo PlO _ P^ 
, (/^oo)^ Poo 


as(/i6) K^Hb) + 0{K^a‘^, fp) .(67) 


V Poo ' V 

The hnal expression for the Wilson coefficients at the low scale is: 

CM= E aPfibTKifibr dr\fib) + 0{alK^) , 

n,m=0 


( 68 ) 


where are functions of only r] = aPuP) / aPjjib) ■, and ratios of the heavy masses. 

At order OialtP), we keep only those contributions to Og and Oio that are proportional to 

The matching conditions, anomalous dimensions and RG-equations presented in Sec¬ 
tions 13.2113.41 do not include the two dipole operators Py^s- For those two operators, it is 
more convenient to consider the so-called effective coefficients 


cfM = c^{pLt) + T.vt\cM)-\cmi^‘0 

1 = 3 '- ^ 

= + 5s(hb) C'f^^''(/i),) + K{^^b) 

-{-as{fib)i^if^b) Gy ^ (pb) -f- 0{ag,K^) , 

Cfi^ib) ^ Cs{^^b) + hi^^b) - \cm] = Gr^''(p, 




(69) 

(70) 


where, in dimensional regularization with fully anticommuting 75 , y = ( 0 , 0 , —— y, — y) 

and 2 ( = (0,0,1, —1,20, — y)- The effective coefficients Gy°°^^^(/ib), ^^^^^(hfe) and 

can be found in Eqs. (10)-(22) of Ref. [33], while Gy°^^®^(/Xfe) can be found in Eq. (12) of Ref. [34]. 

Following Ref. [9], we take into account the complete 0{asK,) term in Cf^{yb) rather than 
only its mi/(M^s^)-enhanced part (as Section [2] would imply). An explicit expression for 
{fib) can be found in Eq. (30) of Ref. [27]. In Tables [3] - [5], we present the relevant 
{fib ). We £x the input parameters to their central values (specihed in sec. |2]) and choose 
fib = [2.5, 5,10] GeV and /xq = 120 GeV. 


4 Matrix elements I 

Once C^"'’^\fib) is found, one needs to calculate the on-shell b —>■ sl~^l matrix elements {Pp 
of the corresponding operators. In the present section, we consider those parts of the matrix 
elements that originate from diagrams with no photons inside loops and/or bremsstrahlung 
photons. These parts are unrelated to the ln(mft/m£)-enhanced correction to the decay width. 
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(00) 

(01) 

(10) 

^(nm) 

[ -0.763 , -0.544 , -0.379 ] 

[-0.180, -0.0835, -0.0378] 

[ 13.764 , 14.943 , 16.066 ] 

^(nm) 

[ 1.054 , 1.029 , 1.015 ] 

[0.248, 0.158, 0.101] 

[ -1.746 , -1.376 , -1.050 ] 

^(nm) 

[ -1.10 , -0.571 , -0.283 ]10-2 

[-1.22, -0.400, -0.125]10-3 

[ 5.28 , 7.98 , 8.38 ]10-2 

^(nm) 

[ -0.113 , -0.0741 , -0.0486 ] 

[-1.62, -0.697, -0.297]10-2 

[ -0.690 , -0.343 , -0.143 ] 

r^(nm) 

'-^5 

[ 1.04 , 0.547,0.274 

[1.17, 0.387, 0.122]10-^ 

[ -1.60 , -1.55 , -1.36 ]10-2 

r<{nm) 

'-^6 

[ 2.32 , 1.17 , 0.563 (lO'^ 

[2.51, 0.801, 0.245]10-^ 

[ -0.656 , -1.92 , -2.17 

^(nm) 

'-^3Q 

0 

[-5.03, -3.72, -2.66]10-2 

0 

^(nm) 

'-^4Q 

0 

[-2.13, -1.04, -0.49]10-2 

0 

'-'5Q 

0 

[-6.08, -1.71, -0.30]10-6 

0 

^(nm) 

^6Q 

0 

[2.12, 1.03, 0.485]10-3 

0 

^{nm) 

0 

0 

0 


Table 3: Numerical values of the relevant {k 7 ^ 7, 8 , 9,10) for = [2.5, 5,10] GeV. 




cr\y^) 

(00) 

(01) 

(10) 

(11) 

[ -0.362 , -0.320 , -0.285 ] 

[ 3.20 , 3.33 , 2.82 ]10-2 
[ 1.625 , 1.171 , 0.690 ] 

[ 4.132 , 4.314 , 4.397 ] 

[ -0.168 , -0.151 , -0.138 ] 


Table 4: Numerical values of the relevant for /i^ = [2.5, 5,10] GeV. 


One-loop penguin contractions of the 4-fermion operators give the following contributions 
to the matrix elements: 


(p^)peng 


M^P,) 


tree 


+ m; 


(Pi) 


tree 


(/^b) 


+ ^/°(-Plo)tree 


(71) 


The above formula holds also for the tree-level matrix element of P 7 , the one-loop matrix element 
of Pq, and for those parts of the two-loop 0{asaem) matrix elements of the 4-quark operators 
where the gluon couples to the closed quark loop. The coefficients Mf are summarized in 
Table [ 6 ] in terms of the functions Ff‘{s) and 


fi{s) = In — + p- (giyc) + ^ In —) + p\g{,yb) + p- (Ins-ivr) + pf , (72) 

Pb V 9 rricJ 

fnS) = 81n^-3(gW + 5ln]][()+5(ln5-in-y. (73) 

Here, i/a = 4m^ pQ[g/m^+£_, the function g{y) is given in Appendix and the numbers p are 
collected in Table [71 The functions Ff'{s) can be found in Eqs. (54)-(56), (71) and (72) of 
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Gr'(h.) 

cir(h.) 

(00) 

0 

0 

(01) 

[ 5.025 , 3.722 , 2.664 jlO-^ 

0 

(10) 

0 

0 

(11) 

[ 2.003 , 1.934 , 1.863 ] 

[ -4.222 , -4.222 , -4.222 ] 

(20) 

0 

0 

(02) 

[ 0.376 , 0.208 , 0.104 jlO-^ 

[ 1.081 , 0.489 , 0.218 ]10-2 

(12) 

[ -6.614 , -4.317 , -2.810 ] 

[ -5.854 , -3.798 , -2.458] 

(21) 

[ 5.645 , 3.538 , 1.193 ] 

[ 5.105 , 6.380 , 7.631] 

(22) 

[ 36.814 , 27.320 , 20.275 ] 

[ -32.014 , -36.090 , -39.764] 


Table 5: Numerical values of the relevant for fib = [2.5, 5,10] GeV. 



M'i 

Ml 

Ml« 

i=l,2 

fi{s) - a^K Ff{s) 

-a'in Fl{s) 

0 

i=3-6,3Q-6Q,b 


0 

0 

i=7 

0 

oisK 

0 

i=8 

P|(s) 

-a^K F^{s) 

0 

i=9 

1 4- asK fr'^is) 

0 

0 

i=10 

0 

0 

1 


Table 6: Coefficients Mf that appear in Eq. flTB . Here, ttg and k are taken at the scale /Xfe. 


Ref. [4] where they are given in terms of an expansion in s up to In the range of s 

that we consider here, the accuracy of these expansions is excellent, as can be seen in Fig. 8 of 
Ref. [6] where the same functions are numerically evaluated for arbitrary s. 

For what concerns the remaining contributions to the NLO and NNLO QCD matrix elements 
of ^ 7 , 9,105 the virtual and real corrections can be effectively taken into account via the following 
redefinitions of the squared tree-level matrix elements in the expression for the decay width: 


I (T9)tree 
I (f^io)tree 
I {P7)tTee 


Re ((P7)tree(R9)tree) 


KP9)treer 

I (Rlo)tree 

I (R7)tree|‘ 


1 8 Cts tl79g^(s) -f- 16 Ctg UJgg 

1+83, wljl„(s)] , 

1 + 8 5 , . 



Re {{P7)tyee{Pil)lree} 1 + 8 tt, 


(74) 

(75) 

(76) 

(77) 


where the functions calculated in Refs. [4,9] are listed in Appendix El 

The remaining contributions to the NNLO matrix elements of the 4-quark operators originate 
from diagrams where the gluon does not couple to the quark loop. Thus, they are given by the 
same functions of s as in Eq. (1741) . 


18 





















Pi 



Pa 

n 


pQ 

pQ 

4 

pQ 

5 

pQ 

-^6 

pb 


4 

3 

1 

6 

0 

60 

0 

4 

0 

40 

0 

0 


0 

0 

7 

2 

2 

3 

-38 

32 

3 

7 

6 

2 

9 

38 

3 

32 

9 

-2 


0 

0 

2 

8 

32 

128 

74 

8 

752 

128 

0 

9 

27 

9 

27 

27 

81 

27 

81 

p* 

16 

4 

2 

8 

136 

320 

358 

8 

1144 

320 

26 

27 

9 

27 

81 

27 

81 

81 

243 

81 

243 

27 


Table 7: Numbers that occur in the four-quark operator matrix elements in Eq. fl72p . 

5 Matrix elements II 

In this Section, we calculate those electromagnetic corrections to the matrix elements of the 
4-fermion operators that are responsible for the ln(mj/m£)-enhanced correction to the decay 
width. In section Em we cover in great detail the calculation of QED corrections to (Eg). In 
Section 15.21 we give the logarithmically enhanced QED corrections to the matrix elements of 
Pi with i ^ 9. 


5.1 Corrections to (Pg) 

Electromagnetic corrections to the matrix element of Pg are infrared divergent and must be 
considered together with the corresponding bremsstrahlung. The dilepton invariant mass dif¬ 
ferential decay width is not an infrared safe object with respect to emission of collinear photons. 
Hence, electromagnetic corrections contain an explicit collinear logarithm ln{ml/mj). The co¬ 
efficient of this logarithm vanishes when integrated over the whole phase space but not if we 
restrict it to the low-s region. 

In this calculation, we adopt the NDR scheme with D = 4 — 2e. The NDR scheme is suitable 
for our calculation since no Levi-Civita tensor survives in divergent terms proportional to 1/e or 
1/e^. Thus, all the Levi-Civita tensors can be evaluated in D = 4 and are therefore well-defined. 

In the hrst step, all the external particles are taken to be on-shell, and, in addition, all 
the hnal state particles are treated as massless (m^ = = 0). This implies that all the 

collinear divergences are dimensionally regularized, and that the collinear logarithm appears as 
a residual 1/e. Later, we will be able to re-express such a residue in terms of \n{ml/mj) using 
the photonic splitting function of the lepton. 

In the next two subsections, we present the calculation of virtual and bremsstrahlung cor¬ 
rections. In the last one, we show how to change the collinear regulator from dimensional to 
the physical mass regularization. 

The calculation involves the following kinematical invariants: Sjj = 2 where i,j G 
{/+,/“, 5,6,7} = {l,2,s,6,g}. 
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Figure 1: Examples of diagrams contributing to the virtual (left) and real (right) electromag¬ 
netic corrections to the matrix element of Pg. 


5.1.1 Virtual corrections 

In order to obtain the virtual corrections, one has to consider one-loop diagrams of the current- 
current type. There are in total six such diagrams, one of which is shown on the left in Figure [H 
The sum of the six amplitudes contains infrared and ultraviolet divergences. The latter cancel 
after the addition of counterterms. The next step is then to compute its interference with 
(F 9 )tree which yields an expression 1 ^ 1 /( 512 , si*, S 2 s). Finally, one has to integrate Ky over the 
phase space. The phase space measure for a three particle massless final state in D dimensions 
is given explicitly in [35]. Since Ky does not depend on angular variables we can immediately 
integrate them out 

dPSs = / dPS^ = Sis, S2s) dsi2 dsis ds2s 

Jn 

2 ” 8+6e^— ^-|-2e 1—2e 

= -ppTr'^yY^(]~—~e) ('^12 Sis S2S) '^(5(1 — Si2 — Sis — S2s) dSi2 dSis (iS2s • (78) 

By means of this expression we obtain the final contribution from virtual corrections via 



5.1.2 Real corrections 

In order to cancel the infrared singularities present in Ty one has to add the corresponding 
bremsstrahlung contribution. There are four diagrams, one of which is shown on the right in 
Figure [H Contrary to the case of gluon bremsstrahlung, the photon couples to all external legs, 
which makes the calculation more involved. The sum of the four amplitudes has to be squared, 
yielding an expression hs, S 2 s, hq, hq, Ssq, Sih), where 


Stri = 1 — Si2 — Sis ~ S2s — Slg + S2q + Sgq 


(80) 
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is the triple invariant. The corresponding phase space measure for the four particle hnal state 
can also be found in [35]. After integration over angular variables, it reads 


dPSi = I dPSi = jd 


6 e 


2-3e 


2 - 12 +lOe ^-5+3e ^^2'^ 

r(|-e)r(l-e)r(i-e) 


dsi 2 dsis dsiq ds 2 s ds 2 q ds 


sq 


X (-A 4 ) 2 ^ ■ 0 (-A 4 ) ■ 5(1 - S 12 - Sis - Siq - hs - hq “ S^g) . (81) 


In the above equation, the Gram determinant is given by 

A 4 = (Si 2 Ssg)^ + (SisS 2 g)^ + (SigS 2 s)^ — 2 (Si 2 SisS 2 gSsg + SisSiqS2sS2q + Si2SiqS2sSsq) ■ (82) 


The phase space measure is completely symmetric in {1, 2, s, g}, but since we stay differential 
in Si 2 we can only make use of the symmetries 1 2 and s gUJ The use of these symmetries 

is, however, essential since the number of distinct terms in gets reduced signihcantly. In 
addition, all terms of the form A/(sigSsg) and A/(s 2 qSsq) as well as B/(siqStri) and B/(s 2 qStri) 
drop out by means of the 1 2 symmetry. 

Another crucial point is to choose for each term in Kj^ the order of integration in a suitable 
way in order to ensure that all terms up to and including order can be found analytically. 
Appendix [C] is devoted to this rather technical issue. The QED bremsstrahlung contribution 
hnally reads 


Tr 


r dPSj 

J dsi2 


Pr(^12} "Sis, S2s} ^Iqi S2qi ^sqi 


(83) 


In the sum of Ty and Tr the 1/e^ terms cancel as well as the part of the 1/e terms, whereas 
the collinear divergences proportional to Qf/e remain. 


5.1.3 Prom NDR to mass regularization 

As we have stated earlier, the differential decay width is not an infrared safe object with respect 
to emission of collinear photons. This means that, as long as the lepton is treated as massless, 
the sum of virtual and real corrections is not free of collinear divergences. If we had kept 
the lepton mass different from zero during the whole calculation, the sum of virtual and real 
corrections would have been hnite. However, the computation of the diagrams and the massive 
phase space integrals in Ty and Tr would have been much more tedious. 

The translation of the 1/e pole into a hi(rnl/mf) corresponds to changing the regularization 
scheme and is complicated by the presence of soft infrared singularities. The correct procedure 
is to start with constructing an observable that is infrared safe and, consequently, regularization 
scheme independent. Only at this point we can switch to the mi regulator and obtain our hnal 
result. As an intermediate step, we construct a differential branching ratio where s is identihed 
as follows: 

^ ^ \(PhPPt2+P-i)‘^l^b Pi W (Ph Ph) (- 34 ) 

\(pii P P(. 2 )‘^/^b otherwise. ^ 

II In the terms containing stri in the denominator, only the 1 2 symmetry can be used. 
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Figure 2: Splitting function kinematics. The photon is emitted by a quasi-real lepton. 


In order to switch to this intermediate observable we must subtract the collinear decay width 
differential in the dilepton invariant mass and add the same quantity but remaining differential 
in the triple invariant. 

The calculation of the differential branching ratio in the collinear limit is done with the help 
of the NDR-scheme splitting function for the massless lepton. The splitting function in this 
scheme can be derived from Refs. [36,37] and reads 


&\x,E) = 4d, 


1 -I- (1 — x)" 


X 


■i + l„^ + l„( 2 - 2 x) 




X 


2x 


X 


, (85) 


where E is the energy of the incoming lepton and xE is the energy of the emitted photon. See 
Fig. [2] for a pictorial view of the kinematics. 

The fully differential decay width in the collinear limit is given by (here and in the following 
we omit the factor 8G‘j^\VtbVts\‘^ stemming from the effective Lagrangian): 


dF^Q{i(si2,sis,S2s,a;) = {2mb) [f!f^{x,Ei) + f!f>{x,E2)\ {{Pg)treel^ dPS3 dx 

= m'i;^ &\x,Ei) \{Pg)tree\^ dPS3 dx , 






( 86 ) 


where a:, S 12 , S 2 s G [0,1], Ei = mb{l — S 2 s )/2 and we used the ii <->■ £2 symmetry. The 
collinear decay width differential in the triple invariant (s = + Pe 2 PP'r)‘^/'^b) reads 


coll,3 -1 / 

ds ~ Jo 


dx dSls dS2sM3{s,Su,S2s) f^"\x,Ei) |(P9)tree 


2 

S12— 


(87) 


The collinear decay width differential in the dilepton invariant mass (s = {pi + XP 2 Y/rnl) reads 


dF 


ds 

where x = 1 


(d 

coll,2 —1 

= mu 


'■ 1 -^ dx 


doc ^ ^ 

— / dsu / dS2sM3{s/x,Sis,S2s) fi\x,Ei) |(P 9 )tree|J^ 2 ^^/s,( 88 ) 

OC ^ 0 0 


X, and the non-linear change of variables S 12 —>■ s/x also implied a distortion of 

the x-integration domain. The addition of —-^ 

removes the remaining e-pole from the differential decay width. 


coll ,3 -results of previous subsections 
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We are now free to convert back this observable to the usual one (in which s is always 
the dilepton invariant mass) using mass regularization. To this extent, we need the splitting 
function in this scheme [36]: 


/f)(x,E)=4«, 


T + (1 — 


X 


E 

In-h ln(2 

mi 


2x) ) —l-\-x- 


x^ 

— In X ■ 
2 


(2 — x)" 


ln(2-x) .(89) 


ip(m) ip(m) 

-^ where dT^'") 

as as 


The original differential decay width is then obtained by by adding 
is obtained from dT^*^^ via —»• Therefore, the total correction term is given by the 

following double difference: 


Ts 

2mb 


dr 


(m) 
coll, 2 


dr 


h) 

coll, 2 


dr 


(m) 
coll,3 


dr 


h) 

coll,3 


ds 


ds 


ds 


ds 


(90) 


Note that only the ^'-independent difference f!f\x, E) — f^\x^ E) enters in the total correction 
term. Hence, we can perform separately the (x, S 12 ) and (si^, S 2 s) integrations. The tree level 
squared matrix element of Pg integrated over the phase space reads 


— 7+6e,^— 2+26 


C^[Sl2 = 


[ml)^ 


ds 


12 S 12 (1 


r(i- 6 ) 

and the total correction term is hnally expressed as 


S 12 ) 


\2-2e 


r(2 — e) 2 . 512(1 — e) -|- 1 
r(2 - 2e) 


3-2e 


r 1 


Tc = 2 


l-s 


Idx 


/P(i) - 


cr s 




X 


X 


(91) 


(92) 


/7 v-y J'y j Q 

LO 0 ^ 

Both integrals in Eq. (j9^ are infrared divergent for x —> 0, but their sum is not. 

The sum Tv + Tii + Ts is now free of divergences and contains an explicit collinear logarithm 
ln(m^/m^). The coefficient of this logarithm vanishes when integrated over .§ 12 . This means 
that if we had considered the total branching ratio instead of the differential one, the sum of 
Tv + Tr would have been already finite and the inclusion of Ts would have become unnecessary. 
However, the coefficient of the collinear logarithm is large and positive for low S 12 and large 
and negative for high S 12 . Furthermore, this term renders by far the major contribution to the 
electromagnetic corrections. In the sum Ty + Tr + Ts the coefficient of Q'^ is up to a color 
factor proportional to the QCD-function Ugg\s) from Eq. (11271) . providing another check for 
our result. Inserting = —1/3 and Qi = —1 finally yields 


'-V 


T 


R 


with 




(em), 
99 I 


Ts = 


= In 


a^ml (1 - S 12 ) (1 + 2si 2) (em) / A \ 

■OOgg (_Sl2j 


ml 


ml 


24 7r3 


1 4s - 8s^ 


(93) 


-Tz2(s) -I--TT^ — 

9 ^ ' 27 


6 (1-s) (l + 2s) 

121 - 27s - 30s^ 


-h ln(l 


72 (1-s) (l + 2s) 


'-3 - 10s - 17s2 + 14 s 3 ^ 17 ln(l - s) 


(1 — 6s^ -I- 4s^) Ins 
^ 2 (1-s)^ (l + 2s) 

(41 -|- 76 s) ln(l — s) 

36(1 +2 s) 

(1 — 6 s^ + 4 s^) In^ s 


18(l-s)^ (l + 2s) 


18 


In s — 


2(l-s)^ (l + 2s) 


(94) 
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The contribution that we have calculated can be effectively taken into account via the following 
substitution: 


\CM{P9) tree | 


|C'9(hb)(-P9)tree|^ 1 + 8 ttg CUgg 


(95) 


5.2 Other log-enhanced corrections 

The QED corrections to the matrix elements of Pi with i ^ 9 contribute to the branching 
ratio at order 0(5^^^). Consequently, following the outline in Section [2l we include those 
contributions that are enhanced by an explicit ln(mb/m£). The relevant terms in the amplitude 
are 


A OC {C 2 + Cf Cl) asH f 2 {s) + Cg (Pg)tree + C'io(Plo) tree T Cy {P 7 ) tree 


(96) 


where the f 2 {S) is defined in Eq. (|72|l . Here we have dropped the NNLO QCD corrections to 
the matrix elements as well as the terms proportional to the small penguin coefficients C'j(Q). 
After squaring and under the assumption that Ci and C 2 are real, we obtain 


|kl| OC ICgl + a^K {C2 + CfCi) |/2('S)| + 2 agK Iie[f2{s) {C 2 + C p Ci)Cg] |(P9)tree 


+ 2 Re C7Cg + asK C'i{C 2 + Cp Ci)/ 2 (s) {P'i)iree{Pg)tree 

+ 1^7^ l(-Pi') tree | TlCior \{Pw) tree | 

The fully differential decay width in the collinear limit now yields 

drS(si 2 ,Si,,S 2 s,a:) = f^^\x,Ei) \A\^ dPSiidx . 


dr 


(m) 


dr 


(97) 

(98) 

(m) 


These corrections are induced by collinear photon emission and are given by —- 

where we retain only the \xi{mi,/me) term in E). The result reads 


dAT 

ds 


Gpuil 


~ ■s)^Q:sK‘\ 8 (1 + 2s) 


+s.KRe (C2 + CfC\)c; + syycj + CfCi)" wS°"(s) 


ic,pa,<r>(s) + ic.„r'=<„';'„'(9) 


2 , ,(em). 


\ 2 (em), 


+96 


a,K Re [CjC;] uj\f^\s) + Re {C 2 + CpCi)C; 


,(em), 


+8 (4 + -)a‘^iA ICyp a; 77 ™^(s) > , 


(99) 


where (Ugg™^ (s) was already found in the previous section. The other cn-functions read: 


^1010 




l + 4s —8s^ , , ... (1 — 6 s^ + 4s^) Ins 

+ ln(l - s) - ^^- - - 

2 (1-s)^ (l + 2s) 


6 (1-s) (l + 2s) 


, ( 100 ) 
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4r'(s) 


LO. 


(em), 
79 I 


(em) 

^29 


(em) 

^^22 


4r’(s) 


In 

In 

In 

In 


rr^ 

nil 


mt 


rrio 


mt 

.ml 


mt 


mt, 


2(l-s) (2 + s) 


+ ln(l — s) — 


§ (—3 -j- 2 3 


2(1 


(2 + s) 


ln(s) 


-^-IT +ln(l 

Si(g)+7S((g) ~ 
8 ( 1 -s) 2(1 + 2 s) 
Som 


(-1 + 2s-2s2) , 

■s) + —tt: -779—- ln(s) 


2(l-5r 


8(1 - s)2(l + 2s) ^ 9(1 - 1)2(1 + 2s) “UgcvJ 


Ei(s) 


In 


f^b \ 


81 


I , .(em)/.x ,2/ ^^b \ 

+ ^^ioio(s) In , 

^3(g) + i ^i(s) 
96(1-s)2 






( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 


The functions Sj have been evaluated numerically in the low-s-region (for fixed values of m^ 
and rUc): 


Si(s) = 23.787- 120.948 5 +365.373 s^- 584.206 (106) 

S{(s) = 1.653+ 6.009 5- 17.080 5^ + 115.880 5^ (107) 

^ 2 ( 5 ) = 11.488 - 36.9875+ 255.330 52 - 812.3885^ +1011.791 5\ (108) 

E3(5) = 109.311 - 846.039 5 + 2890.115 5^-4179.072 5^ (109) 

S^(5) = 4.606+ 17.650 5- 53.244 5^ + 348.069 5^ (110) 


6 Collinear logarithms and angular cuts 

The explicit logarithm of the lepton mass signals the presence of a collinear singularity whose 
appearance in the differential branching ratio is strictly related to the definition of the dilepton 
invariant mass. As explained in Sec. 15.1.31 this logarithm disappears if all photons emitted by 
the final state on-shell leptons are included in the definition of s: (p^j (Vh +P ^2 +P 7 )^- 

Let US consider a cone (of angular opening 9) around an on-shell lepton with momentum 
Pi. For all photons emitted in this cone we have: < (pi + p^)^ < ~ 2i7|(l — cos 6*), 

where Ei is the energy of the lepton (usually of order mb). Consequently the effect of including 
such photons in the reconstruction of the lepton momentum can be roughly approximated by 
replacing by some scale of order A in the collinear logarithm. 

Both Babar and Belle include sufficiently collinear photons in the definition of the lepton 
momentum. However, the imposed cones are so narrow that the effect for the muons is negli¬ 
gible, i.e. the separation of muons and collinear photons is practically perfect [10]. Thus, our 
expressions containing ln(m^/m^) are directly applicable in this case. 
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For electrons, the situation is more complicated. In both experiments, the cone is dehned in 
the laboratory frame and has polar and azimuthal angles around 45 mrad and 5 mrad, respec¬ 
tively. Hence, A is of the same order as m^, which makes the QED corrections for the electrons 
similar to those for the muons. Nothing more precise can be said without applying dedicated 
Monte Carlo routines that would take into account the experimental setups in detail. 


7 Formulae for the branching ratio 

In Section O we have expressed the branching ratio in terms of the quantity ^u{s)/^u- In 
the present Section, we express this quantity in terms of the low-scale Wilson coefficients and 
various functions of s that arise from the matrix elements. The main formula reads 




$ 


u 


i<j 




( 111 ) 


where ^ only for i = 7,8 (see Eqs. flUUj) and fl7n|) b The functions s) can 

be expressed analytically in terms of the coefficients Mf listed in Table [6] and of the following 
building blocks 


S99 = (1 - s)^(l-h 2s) |l8 ttg a;99^(s)-F 8 OsK a;99“^(s) 


+ 16 al (s) + + 4 m(^)49^(s)] } 

+ 6 ^ (1 - 6s2 + 4s3) , 
mi 

S77 = (1 — s)^(4 + ^) |l + 8 dg a;77^(s) + + 8 dgK a;77™^(s)| 

+ 24 ^ (2s2 - 3) , 
mi 

S79 = 12(1 — s)^ |l + 8 5g a;79^(s) + + 8 dg/t a;79™^(s)| 


+ 24 ^ (1 - 6s + 4 s2) , 
mi 


Sioio — 5*99 + 8 Q!sfi:(l — s)^(l + 2s) ct;|oio (s) — c^gg '’(s) 


,(em) / 


( 112 ) 


(113) 


(114) 

(115) 


The functions are listed in Appendix |Al The functions have been given in Eqs. (|M|1 
and (I100l) - (ll05p . The numbers = (dvr^ —25)/12 and ~ 27.1+/5 om^^^ ln(/if,/mb) originate 
from the QCD corrections to 6 —*• Xuev [38], while the quantity ^ (j^-i _ 1 ) stands 

for the logarithmically-enhanced QED correction to this decay [39]. The Saa include non- 
perturbative Cl(l/m^) corrections that one hnds by taking the limit m^ —> 0 in Eq. (18) of 
Ref. [22]. 
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The explicit expressions for the functions Hij read 


E |M^|2 + Re{MjMf*) Sj, + Ai/,, , when i = j 


H^,= 


E 2MfMf* Saa + [mJM]* + MfM]*) Sjg + AHij , when i ^ j 


(116) 


[^A=7,9,10 

It is assumed that all the products in Eq. (11111) are expanded in 5s, n and A 2 , and that higher 
orders are neglected (see Section [2]). The quantities 


AH{j hij T Cij T Cij 


( 117 ) 


that need to be included only for z = 1,2 stand for additional bremsstrahlung ( 6 p), non- 
perturbative 0{l/ml) corrections (qj) and additional ln(mf/m£)-enhanced electromagnetic 
corrections (cij). Specifically, the non-vanishing Cij that we include read 

622 = 8 (1 - s)^ (12s) 5 sK^ a;22“^(s) 

627 = 96 (1 — s)^ 5g/7^ a;27™^(s) 

629 = 8 (1 - s)^ (12s) 5g/7^ a;29“^(s) 

Cll = y 622 

612 = I 622 

eij = | 62 j, forj = 7,9. (118) 


The 0{l/m^) non-perturbative contributions were calculated in Ref. [21] 


C 27 

~2 2 *^^2 / 

C 29 


C 22 

= “-NmA 

Cll 

= C 22 

C 12 

= 1 C 22 

Cij 


where r = 

1/2/c = 


, 1 -f 6s — s^ 


ReF(r) , 


(1 — s)^(2 -I- s) ReF(r) , 


for j = 7, 9 


(119) 


The finite bremsstrahlung contributions bij appear at NNLO in Ref. [5], where the notation 
is very similar to the one proposed here. We do not present these corrections here but do 
include them in the numerical analysis. 


27 








8 Acknowledgments 

We would like to thank Piotr Chankowski for bringing the problem of electromagnetic loga¬ 
rithms ln(M|f/M|) to our attention. We are grateful to the authors of Ref. [9] for a careful 
reading of the manuscript and useful remarks. We are deeply indebted to Thomas Gehrmann 
and Aude Gehrmann-De Ridder for many discussions and help on the phase space integration. 
We thank Jeffrey Berryhill and Akimasa Ishikawa for information concerning the treatment 
of collinear photons at BaBar and Belle. This work was supported in part by the Schwei- 
zerischer Nationalfonds. M.M. acknowledges support from the Polish Gommittee for Scientihc 
Research under the grant 2 P03B 078 26 and by the European Gommunity’s Human Potential 
Programme under the contract HPRN-GT-2002-00311, EURIDIGE. Research partly supported 
by the Department of Energy under Grant DE-AG02-76GH03000. Fermilab is operated by 
Universities Research Association Inc., under contract with the U.S. Department of Energy. 


A Various functions 

The loop functions that appear in the text are: 
A{x) -- 


—3x^ + 2x‘^ , 8a;^ -|- 5a;^ — 7x 

mx 


Y{x) 

W{x) 

S{x) 

Xix) 

E{x) 


2{x-l)^ 
3x'^ 


Inx + 


12(a;- 1)3 
x"^ — 4:X 


8 (a;-l) 2 ‘“‘^ ' 8 (a;-l)’ 
-32a;^ + 38a;3 + - 18a; 


18(a;- 1)4 


In a; 


-18a;4 + 163a;3 - 259a;2 + 108a; 

36(a;- 1)3 ’ 


3a;4 


2 (a;- 1)3 

3a;^ — 6 a; 
8 (a;- 1)2 


Ina; + 


In a; 


X 


llx^+ Ax 


A{x — 1)2 


x^ -I- 2 a; 
8 (a;- 1 )’ 


a;(18 — 11a; — a;^) a;^(15 — 16a;-|-4a;^) 


+ 


Ina;-Ina; 

3 


12(1-a;)3 6(l-a;)4 

The following function appears in the matrix elements of the 4-quark operators: 

— in, when y < I, 


giy) = - ^{2 + y)^\l - y\ 


1 

j In 

i 

I+a/T+J 

I 1 - 1/11 

l-Vl-y 

( 2 arctan 


y-l ’ 


when y > 1 , 


( 120 ) 

( 121 ) 

( 122 ) 

(123) 

(124) 

(125) 

(126) 


The functions that include the sum of infrared divergent virtual and real contributions 
to the matrix elements of P 7 , Pg and Pio are: 


u, 


99 ' 


4 2 

s) = --L? 2 (-s) - - ln(l - s) Ins 


5 -I- 4s 


:n 


■ ln(l 


3(1 +2s) 




2s(l + s)(l-2s) ^ 5 + 9 s-6s2 

3(1 - s)2(l + 2s) ^ 6(1 - s)(l + 2s) ’ 

(127) 


= -19.2 + 6.1s+ (37.9 +17.2lns)s2-18.75^ , 

(128) 

‘^1010 ('5) 

= iOgg (s) , 

(129) 

<^77^ (s) 

- +' 3(2 + .-)'“*^ 



2s(2-2s-s2), ^ 16-11s-17s 2 8 , ,/it, 
3(l-s) (2 + s) 18(l-s)(2 + s) 3 ^rrib 

(130) 

^7% (s) 

= ^Li 2 {s) ^ln(l s)lns +7^^n(l s) 



2 s (3 - 2s) 5 - 9s 4 /Xfe 

- 15 - In s H-;- -7 -In — . 

9(l-s)' 18(l-s) 3 

(131) 


The function Ugl\s) has been extracted [25,26] from the 0(as) corrections [40] to the semilep- 
tonic decay. The functions a;j 7 ^(s) and have been calculated in Ref. [4]. Note that 

a; 77 ^(s) in the s —> 0 limit reproduces the 0(as) correction [41] to the matrix element of Py in 
the b Xs'y decay. The function ujgg\s) was extracted [9] from the corrections [42,43] 

to the spectrum of the b —> X^ev decay. The approximate formula in Eq. (I128p is valid in the 
range 0 < s < 0.4. 

The function F{r) that arises in the the non-perturbative corrections reads [21] 


F{r) 


3 

2 r 


, arctan - 

V^r(l-r) V 


2 ^ r(r-l) \ l+Y^l-l/r 


- 1 , 

+ m 


when 0 < r < 1, 
1 , when r > 1. 


(132) 


B V and 


The numerical diagonalization of the matrix IRdo) yields: 


OjI 


- 1.04348, -0.899395, -0.521739, -0.521739, -0.422989, 0.408619, 

0.26087, 0.26087, 0.26087, 0.145649, 0.130435, 0,0 


(133) 
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and 


V 


0 

0 

0.942522 

0.0253179 

0 

0 



0 

0 
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-0.0084393 

0 

0 
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0 

0 / 
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C Details of the bremsstrahlung calculation 


This last appendix is devoted to some technical details of the bremsstrahlung calculation. We 
will integrate the sample kernel Kr = over the four particle phase space, show how the 

Gram determinant factorizes, and explain how to extract all terms up to and including order 
e° analytically. Omitting bothersome prefactors and, in addition, removing all hats from the 
invariants Sij we consider the expression 
1 

A := dsi2 j dSisdSiqdS2qdS2sdSsq5{l — Si2 — Sis — Siq — S2s — S2q — Sgq) 

0 

X (-A4)-^-^ 0(-A4)W^2/ , (135) 

where the Gram determinant is given by 

A4 = (Si2<Ssq)^ + + (SigS2s)^ — 2 (Si2SisS2gSsg + SisSiqS2sS2q + ’Sl2<SigS2s<Ssg). (136) 


30 



The first integration is over the ^-function. It is done by means of the variable Sgq, yielding 


A 


1 — S12 1 — Sl2 — Sls 1 —■512—«ls —Slq 1—512 —■51s—■5lq —.529 



0 0 0 0 


X 


(-A 4 ) 


^-^0(-A4) 


Ssq = l —512—Sis—5I9—S29 —S2s 


(137) 


Substituting .s^g = 1 — ,§12 — Si^, — .Sig — S 2 g — S 2 s in the Gram determinant yields an object that 
can be transformed into a quadratic polynomial in either of the variables sig, S 2 s or S 2 g, he. 
in either variable that does not accompany Sgg in the quadratic piece of the Gram determinant. 
We choose this quadratic polynomial to be in 52 ^: 


- A 4 


~(Sl2 + 


S 2 s -\- ‘2 B S2s + c 


(Si2 + — S2s){s2s — S2s) , 


(138) 


where sfg are the roots of the quadratic polynomial: 

4. = -B± VB^-C = -B±Ve. (139) 

The ©-function now requires these roots to be reaQ which is equivalent to the condition S > 0. 
From the latter inequality we conclude 


S2q ^ Z {1 — S 12 — Sis ~ Sig) 


with 


S12 + Slg 
S12 + Slq -|- Sis 


< 1 


(140) 


The above roots now fulhll the inequality 0 < S 2 s < s^s A 1 — S 12 — sis — sig — S 2 g which leads 
to the following new limits of integration: 


1 — 512 1 — 512 —51s 


A = 


\-l-2e 


dSi2 J dSis j dSigs4(Si2 + Sig) 

0 0 

2(l-5l2-5ls-5l9) 5 + 

X J dS 2 gS 2 g J dS2s (sj, - S2s)~^~" (S2s “ Sss)"^"" 


(141) 


Substituting S 2 s = (s 2 s~S 2 s)x + S 2 s, the subsequent y-integration can be done trivially in terms 
of T-functions. 

As a general strategy for the choice of the order of integration, we suggest the following. 
The variable of the hrst integration (^-function) must not be contained in the term of iF/j 
that one considers. If possible, this term should also be free of the variable that one uses to 
factorize the Gram determinant (s 2 s in our case). If the latter is not possible as for instance in 
= SisS 2 sSig^S 2 g^, one should at least factorize the Gram determinant in a variable that does 
not appear in the denominator of This procedure ensures that the hrst two integrations 
can be done in terms of T-functions, and it avoids hypergeo metric functions to emerge at this 
stage of the calculation. 

**otherwise —A4 is negative for all S2s 
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The choice of the subsequent order of integration is governed by the aim to extract all 
divergences as early as possible, this being the reason why we solved the condition S > 0 for 
S 2 q. We now substitute S 2 q = z(l — S 12 — si* — sig) t and perform the t-integration, yielding 
again only T-functions. After simplihcation, we obtain 

1 —SI2 1 —S12—Sis 

A= - dsi2 S 12 J dsis Sig J dSiqSig '^(<§12 + Siq) ^ (<§12 + Sig + — S 12 — — Sig) ^*^.(142) 

0 0 

We now proceed as follows 

111 

• Perform an expansion into partial fractions via — - -- =- - -- 

'Slg('Sl2 + -Siq) <Si2<Sig <§12 (<§12 + -Slq) 

• Snbstitnte = (1 — S 12 — <Sis)(l — u). The n-integration can be carried ont in the first 
term of the above expansion. 

• Snbstitnte = (1 — Si 2 )(l — v), again the n-integration can be done in the first term. 


One now obtains the following expression: 


.4 = 


i r(2-4e) 


2 Fi{—€, 1 — 2e; 2 — 4e; 1 — S12) 


1 1 


— (1 — S 12 ) J du J dv 


u — u) ^"^(l — v) 


[1 - (1 - si2)uvY 


0 0 


S 12 + (1 - (1 - U) 

We now carry out a two-dimensional variable transformation from [u, v) to {y, w) via 


(143) 


y = 1 — uv and yw = v {1 — u) . 


(144) 


The w-integration can now be performed, resulting in another hypergeometric fnnction. After 
nsing the Knmmer relation, also the ^-integration can be done. The hnal result for A reads 


A = - dsi2S 


-l-e 


12 *12 


(1 


- 2 n(-e. 1 - 2 £; 2 - 46; 1 - 


1 


r(2 - 4e) 


S 12 ) 


, , r2(l-e)r(l-2e) .61 

■(1 — S12)--2-^i(l — e, 2 — 2e; 3 — 4e; 1 — Si2)| • 


(145) 


All the divergences have now been extracted in terms of poles and T-fnnctions. The remaining 
task is now to expand the hypergeometric fnnctions in e. This can can be done by means of 
the Mathematica package HypExp [44]. 
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